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Abstract 


Max, Timothy A.; Schreuder, Hans T.; Hazard, John W.; Oswald, Daniel D.; 
Teply, John; Alegria, Jim. 1996. The Pacific Northwest Region vegetation 
and inventory monitoring system. Res. Pap. PNW-RP-493. Portland, OR: U.S. 
Department of Agriculture, Forest Service, Pacific Northwest Research Station. 
PP {0}, 


A grid sampling strategy was adopted for broad-scale inventory and monitoring of 
forest and range vegetation on National Forest System lands in the Pacific North- 
west Region, USDA Forest Service. This paper documents the technical details of 
the adopted design and discusses alternative sampling designs that were considered. 
A less technical description of the selected design will be given elsewhere. The grid 
consists of a regular, square spacing with 5.47 kilometers (3.4 mi) between grid 
points. The primary sampling unit (PSU), established at each grid sampling point, 
consists of a circular, 1-hectare (2.47-acre) plot. The PSU is subsampled with a 

set of different-sized fixed-area subplots, as well as line transects, to assess all 
components of vegetation. The design is flexible and can be used with many 

types of maps. The theory of point and change estimation is described, as well 

as estimates of variation that assess the statistical precision of estimates. 


Keywords: Sampling, plot design, fixed-area plots, line intersect sampling, 
monitoring, National Forest System, Pacific Northwest. 
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Introduction 


In the past, inventories of National Forest System lands (NFS) have been designed 
primarily for estimating potential outputs of and managing the forests for wood 
products. Some information on other resources was collected, but it was often in- 
sufficient to answer complex questions about these resources. In addition, profes- 
sional groups other than timber managers have conducted surveys of other resources 
such as wildlife habitat, soils, recreational features, water, and so forth. These efforts 
were rarely integrated into a comprehensive resource inventory system, plus inventory 
and monitoring were rarely considered together. 


National Forest System lands are managed for a multitude of objectives including 
providing a wide variety of commodities, uses, and scenic and social values that 
people derive from them. Inventory and monitoring systems must be designed to 
provide comprehensive information about the land and vegetation to answer a 
broad spectrum of existing and newly developing questions. For example, infor- 
mation needs on the occurrence, abundance, and distribution of Pacific yew (Taxus 
brevifolia Nutt.) or the occurrence of habitat suitable for the northern spotted owl 
(Strix occidentalis Caurina) would not have been anticipated 10 years ago. But 
such information is much more likely to be available in a broad land and vegetation 
database. Such broad inventory and monitoring schemes also must meet the 
informational needs of the USDA Forest Service (FS) for periodic national and 
Regional assessments, Regional and national planning, and developing and 
implementing forest plans for NFS lands. 


In October 1992, the Washington Office of the FS ordered its western Regions 

to establish a grid of plots across all NFS lands that would be compatible with the 
Forest Inventory and Analysis (FIA) grids. In April 1993, the FS established policy 
implementing ecosystem management of NFS lands. In June 1993, the Washington 
Office directed the FIA units to implement mapped plot designs for all FIA plots in 
the United States (Robertson 1993). The key idea of this direction was to locate 
plots exactly on the ground as planned and map them for conditions of interest. 
Moving some subplots into the same condition class as that of the central sub- 

plot is not permitted. This direction was expected to be fully implemented in 1994. 


Anticipating a need for a scientifically defensible yet practical design for ecological 
purposes, the Pacific Northwest Region of the NFS (Region 6) convened a panel 
of survey sampling experts who met in a series of five meetings to plan such an 
inventory and monitoring (IM) system. This paper summarizes the deliberations 
and conclusions of this advisory group, which recommended a comprehensive IM 
system for Region 6. We believe this is a key step in inventorying and monitoring 
forest land for ecological purposes. Under current monitoring, only large-scale 
change is estimated by using sampling units that are periodically remeasured. 
Various alternatives for implementing remeasurement, such as an annual forest 
inventory system (AFIS), ’ will be explored later. Professional and disciplinary 
groups, whether monitoring watersheds or other areas of interest, can use the 
database developed here as a starting point for their more specific monitoring 
projects. 


Y Hansen, M.H.; Hahn, J.T.; Schreuder, H.T.; Befort, W.A. 
AFIS-1: implementation of an annual forest inventory system. 
Manuscript in preparation. 


Purpose and 
Objectives 


Developing Forest Plans 


Assessing Resources 


Special Studies 


The purpose of this paper is to present a broad-scale vegetation inventory and mon- 
itoring system designed to meet the present and future resource information needs of 
Region 6 and the FS. The system must be integrated by capturing basic land infor- 
mation and describing vegetative structure and productivity in an ecological frame- 
work. In addition, the system must be responsive to changing needs for spatial and 
temporal information at various scales and levels of resolution. 


The objectives of the system are to develop and maintain a defensible database of 
ground resource estimates that can be used to periodically assess vegetation and 
monitor change in vegetation over time. The data are compatible with any map base 
developed for NFS lands in Region 6. The database will be used to develop forest 
plans and assess resources, special studies, and impact analysis. Each objective is 
described in more detail below. 


For monitoring, it is important to note that we have used the term “monitoring” to 
mean assessing broad-scale changes in resources over time. The scope is strictly 
broad-scale; that is, applying to large areas. Many specific types of monitoring are 
defined and used by the FS and other agencies (for example, implementation 
monitoring and effectiveness monitoring), and they are not part of the scope of 

this paper. A good discussion of inventory and monitoring, as well as definitions of 
specific types of monitoring, is given in MacDonald and others (1991:6-8). In this 
paper, we use the term “sampling strategy” to define sample design plus estimation, 
where sample design refers to the approach used to select the actual sampling unit 
locations (Schreuder and others 1993). 


The panel had one critical objective, which was to develop an inventory and mon- 
itoring system that would meet the needs of forest planning and forest management 
decisionmaking. Major components had to include methods for acquiring repre- 
sentative data, methods for summarizing data to obtain estimates, and associated 
database support. Planners now need estimates of both the current status of 
resources and changes in resources over time to make planning decisions. The 

new process will be based on incorporating spatial information for current vegetation, 
potential vegetation, land form, and soils at both Forest and multi-Forest levels. 


Another key objective of the system was broad assessment of resources. Such 
assessment provides consistent resource estimates for addressing subregional, 
regional, and national issues (for example, information required by the national 
Resources Planning Act [RPA] 1974) and an initial source of information to Ranger 
Districts for addressing management issues. Such issues normally are driven by 
users’ concerns for the land and for its associated vegetation, both current status 
and change over time. 


The system had to be flexible to allow for implementing broad-scale special studies 
compatible with the overall IM design. Examples of such special studies are deter- 
mining the pattern of root rot occurrence, the frequency and extent of insect epi- 
demics, and the quantity and spatial distribution of live vegetation and woody debris 
as it relates to forest succession. Although such special studies are observational in 
nature, they provide an effective tool for addressing some issues. 


Impact Analysis 


Review of Literature 


Sampling Designs in Use 


An objective becoming increasingly important is impact analysis. Such analyses 
investigate potential relations among response variables, management practices, 
environmental factors, and other anthropogenic stresses. The acquired vegetation 
databases, both spatial and temporal, must be useful for such analyses. Evaluations 
may include documenting effects of catastrophic events, determining effects of land 
form on management practices, and monitoring change in vegetation and under- 
standing why such change occurs. 


Surveys traditionally have been planned to provide estimates of current status. 
Such estimates have been so useful to management that additional objectives have 
evolved. The original objective of FIA inventories, for example, was to estimate the 
area of forest land by type, stand size, ownership, site quality, and stocking and 

to estimate merchantable wood volume by tree species and diameter class. Over 
time, additional information related to change over time was added, such as growth, 
mortality, timber removals, and success of regeneration. 


Recently, interest has expanded in relations established by using natural resource 
survey data. Establishing such relations and drawing conclusions based on them is 
termed “analytical inference.” Such inferential use of sample survey data requires 
reliable data. Conclusions drawn by scientists from a comparison of the 1982 and 
earlier FIA surveys of Georgia and Alabama is a striking example of such inference. 
In that study, average growth rates in tree diameters decreased (Bechtold and others 
1991, Ouyang and others 1993, Ruark and others 1991) for a screened data set, 
but the cause is still not known. Various hypotheses have been proposed, including 
regional changes in average stand structure of forests, anthropogenic stresses, and 
weather. As a result of the ensuing controversy, efforts were made to determine 
whether survey data could suggest cause-effect relations, and how they might 

more effectively serve analytical uses. 


In the following subsections we briefly review several topics: sampling designs in 
use, estimation and analysis, analytical uses of survey data, and quality assurance 
of data collected. 


Around the world, many broad-scale forest inventories are based on a systematic 
arrangement of grids (Kohl 1994). These designs are adopted primarily for practical 
reasons of implementing and revisiting plots well distributed over a vast area. Three 
well-known examples in the United States are the FIA program, Forest Health Mon- 
itoring program of the FS, and the National Wetlands Inventory of the U.S. Fish and 
Wildlife Service. 


All FIA units use double sampling for stratification. A large number of grid locations, 
either single points or a cluster of points, is located, typically on available photo- 
graphy, and locations are classified into strata (the simplest being forest and non- 
forest) (Birdsey and Schreuder 1992). Permanent, unobtrusively marked sample 
plots are then located on the ground on a grid, typically consisting of a cluster of 5 or 
10 variable-radius plots sampling a 0.4-hectare (1-acre) plot. These plots are then 
measured for variables of interest, which are sometimes multiresource in nature. At 
present the timber variables receive the most attention (Birdsey and Schreuder 1992; 
Hahn and others, 1995). The Forest Health Monitoring program of the FS and the 
Environmental Protection Agency (EPA) uses a cluster of four fixed-area circular 
subplots to sample a circular 1-hectare plot; 0.0168-hectare subplots (radius = 7.32 
meters or 24 ft) are mapped for forest conditions (Messer and others 1991; Scott 
and Bechtold 1995). 


Estimation and Analysis 


Since 1974, the U.S. Fish and Wildlife Service has inventoried the Nation’s wetlands 
through its National Wetlands Inventory project. The statistical design implemented in 
the lower 48 States can be used to obtain reliable estimates for individual states or 
other geographical areas by intensifying the grid (Hall and others 1994). 


Designs are given in Schreuder and others (1993) for sampling recreational use, 
wildlife populations, regeneration, understory vegetation, and timber variables. More 
practical, objective methods are needed for obtaining information on vegetation pro- 
files (Schreuder and others 1993) and the biomass of range vegetation (National 
Research Council 1994). 


Traditionally, interest in NFS lands has centered on the area and condition of the 
resources and how these change over time. Scott and Bechtold (1995) and Williams 
and Schreuder (1995) present the relevant estimation theory for mapped designs. 
The notation used is: 


A = Total area in population (hectare). 

F = Number of condition classes of interest. 

t = Time period t (t = 1, 2). 

Att = Area in condition class f at time t (f = 1, ..., F; t= 1, 2) (hectare). 

Yt = Total for variable of interest at time t. 

AY =Total change in variable Y from time 2 to time 1. 

Ytf = Total for variable of interest in condition class f at time t (f= 1, .... F; t= 1, 2). 
k = Number of subplots in a cluster (k = 5 or 10 depending on FIA unit). 

M = Number of plots in population. 

m = Number of sample plots (clusters). 


BAF = Basal area factor in variable radius plot (VRP) sampling = square 
meter/hectare (or ft7/acre) of basal area represented by a tree when tallied ata 
single subplot. 


Ilti = Probability of selecting tree i at time t (t = 1, 2). 

rs = Radius of small fixed-area plot. 

tr, = Radius of large fixed-area plot. 

Yti = Value of variable of interest for tree i at time t (t = 1, 2). 


vin = Estimated total for variable of interest y on sample plot (cluster) h (h = 1, ..., 
m) at time t (t = 1, 2). 


aj(ft) = Approximate area in condition f (= a, if only one condition occurs) at subplot j 
and time t (t=1, 2). 


dti = Diameter at breast height of tree i at time t. 
gti = Basal area of tree i at time t. 
i = Used to sum over trees on each subplot. 


j = Used to sum over subplots. 


h = Used to sum over clusters of plots. 
as = Area of small fixed-area plot (hectare). 
a, = Area of large fixed-area plot (hectare). 


Estimating area and change—Estimators giving areas by condition classes [A(f)] in 
the population and change in these areas over time [AA(f)] are: 
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Estimating amount and change—The basic estimator of the total amount Y of any 
resource at a given time is: 


Ym / T= ¥n /m (7) 
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For estimating change, the following formulae are recommended. 


Overall change AY — 
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where AY,(f),Y5;,Y2 and Y,are the estimated total change for variable y in condition 
f, estimated totals at time 2 and time 1 in variable y in condition f, and estimated 
totals for time 2 and time 1 for variable y for plants not in the sample plot at time 

1 but present at time 2 (ingrowth, |), respectively. 
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where i(f), S,(f), M(f), and A(f) are estimators of ingrowth, survivor growth, 
mortality, and removals respectively for condition f. If additivity of the estimates 
is not required, a more efficient estimator of overall change A Y is: 


AY2 =AS,-C-M+l, (14) 
where 
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Analytical Uses of 
Survey Data 


The classical variance estimator for these parameter estimates is (Williams and 
Schreuder, 1995): 
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Vi=S(in-¥) /[n(m—1)]. (16) 


h=1 


In addition, Williams and Schreuder (1995) give formulas for estimating the percent- 
age of plots with more than one condition class, the length of the boundaries between 
the condition classes, corresponding changes in these parameters over time, and the 
relevant variances. Instead of using the classical variance estimators as described, 
an alternative is to generate a large number of bootstrap estimates, which allows 
more reliable confidence intervals to be constructed (see Ouyang and others 1992). 


Skinner and others (1989) consider P. Lazarsfeld as the key figure influencing modern 
developments in the analysis of surveys. Lazarsfeld viewed survey analysis as being 
primarily concerned with relations between variables. He was primarily interested in 
methods for causally analyzing contingency tables. 


Deming (1975) distinguished between enumerative (descriptive) and analytical studies 
as follows: In an enumerative study, the number of units of a population belonging to 
a certain class is estimated, whereas in an analytical study, a basis for action on 

the causal system or process to improve the result or effect has to be found. In 

other words, interest in enumerative surveys is in estimating population parameters, 
whereas in analytical surveys it is in explanatory (or predictive) models (Skinner and 
others 1989). Analytical uses of survey data then involves model building and model 
testing. 


To further clarify the distinction between enumerative and analytical studies, consider 
the case where a 100-percent sample of a population is collected. In an enumerative 
study, such a sample provides the complete answer to the question posed, subject 
to the limitations of the method of investigation; that is, an enumerative survey has 
the property that if all N units are observed without error, population parameters are 
determined exactly (Smith 1983a, 1983b). Such populations are called finite popu- 
ations by statisticians. In contrast, inference on the effect of a treatment on a 100- 
percent sample of a stand of trees is still inconclusive in an analytical problem 
because this stand represents an often ill-defined, essentially infinite population 

of such stands. Such populations are called superpopulations by statisticians. 


Traditional sample surveys have been used primarily for enumerative purposes; 
that is, estimation for finite populations. But such surveys are often also used for 
analytical purposes by adopting a superpopulation model; this is vividly illustrated 
by Brewer and Mellor (1973) and Holt and Smith (1976). The alternative uses of 
survey data emphasize the need to recognize that it is the inference to be drawn 
that is inherently enumerative or analytic, and not the survey itself. 


In some of the social and demographic sciences, the enumerative uses have been 
subordinate to the analytical ones (Kendall and Lazarsfeld 1950). In forestry and 
other renewable natural resource disciplines, the reverse is usually true. 


Analysis of data from large, complex sample surveys can be quite difficult (Landis 
and others 1982). Difficulties result from complications introduced by unequal prob- 
ability sampling, imputation for nonresponse, regression estimation and use of 
several levels of sampling and poststratification, and conflict with the assumptions 
made in classical statistical analysis of mutually independent sample observations 
drawn by simple random sampling. In addition, survey data may contain substantial 
and confounded measurement errors, although classical models assume no or 
simply distributed observational errors (Simmons and Bean 1969). There also is 
the question of what superpopulation the sample represents, which is even more 
pertinent in classical statistical experimentation. 


Beyond developing explanatory or predictive models, scientists and managers are 
also keenly interested in cause-effect relations. In considering the conditions nec- 
essary to establish a cause-effect relation, and how difficult is it to achieve these 
conditions, Mosteller and Tukey (1977) note that two of three criteria (consistency, 
responsiveness, and mechanism) must be satisfied. Consistency implies that the 
presence and magnitude of the effect (y) is always associated with a minimal level 
of the suspected causal agent (x). Responsiveness is established by experimental 
exposures to the suspected causal agent and reproducing the symptoms. Mechanism, 
as it relates to the subject of this paper, is established by demonstrating a biological 
or ecological process that causes the observed effect. Only consistency can be 
confirmed by observation alone. Hence, surveys are useful primarily in identifying 
potential cause-effect relations. 


There are many examples, particularly in epidemiology, of falsely claiming to have 
established cause-effect relations. Epidemiologists sometimes cannot conduct 
experiments on people and cannot randomly assign healthy people to exposure 

to potentially noxious substances. Because of these limitations, medical research 
often has not followed scientific principles. Recently, however, rigorous clinical trial 
methods have enjoyed a major resurgence in medical research. Feinstein (1988) 
advocates the following principles for observational studies in epidemiology: stipulate 
a research hypothesis prior to analysis, study a well-specified cohort having a stat- 
istical factor in common (for example, a tree diameter-class), collect high-quality data, 
study possible explanations, and avoid personal biases in detecting possible relations. 


Although inventories in forestry traditionally have emphasized point and change 
estimation, analytical uses of inventory data are becoming more important. An 
example is the effect of treatments on a specific resource. This is particularly likely 
to be of importance for NFS inventories where there is a need to examine the effect 
of treatments applied under operational conditions. Note that effective analytical use 
requires reliable point and change estimates. 


Quality Assurance 


Use of survey data either to identify potential cause-effect relations or to establish 
them is controversial (Schreuder and Thomas 1992). For sample plots, there are 
still many uncontrollable variables (weather, insect epidemics) and essentially 
unmeasurable variables (rainfall, pollution), so many explanations are readily 
available for an observed effect (Such as a change in the population of a bird 
species). Nevertheless, there is no alternative but to use survey data to identify 
possible cause-effect relations. Schreuder and McClure (1991) suggest modifica- 
tions to FIA designs to improve detection of change and identify its possible causes. 
The broad objectives underlining the modifications to FIA proposed by Schreuder 
and McClure (1991) are to generate descriptive statistics, detect changes in such 
statistics, and analyze the data to identify possible cause-effect relations. But 

even with these modifications, the emphasis in FIA surveys would still be to collect 
descriptive statistics efficiently, especially for timber. The fact is that it is difficult 

to even identify reasonable potential cause-effect relations, let alone establish 
cause-effect in the natural resource fields. Simplicity in survey design is more 
likely to lead to success in this regard (Schreuder and others 1993). 


Quality assurance refers to the quality control imposed on the collection and proc- 
essing of data to obtain the most reliable data set and statistics possible for the 
amount of money spent on the survey. The most detailed quality assurance program 
for natural resource inventory and monitoring activities probably is the one used by 
Forest Health Monitoring as required by the Environmental Monitoring and Assess- 
ment Program (EMAP) of the EPA (U.S. EPA 1992). Forest Health Monitoring 
embraces the concept of total quality management (TQM), which is a process of 
continuous improvement and innovation led by managers. Total quality management 
fully integrates management philosophy, planning, and operating methodology. The 
philosophy behind TQM is aimed at achieving total employee commitment to quality. 
To achieve this, TQM focuses on: 


e Identification of customers of the survey. All customers must be identified, 
and effective and continuous communication of their requirements must be 
maintained. 


e Standards and performance. Proactive rather than reactive measures of 
performance must be used. 


e Leadership commitment. This commitment is essential and can be documented 
well by establishing a TQM “culture” during training. 


e Employee recognition. A key ingredient to the success of TQM is to establish 
and implement criteria and mechanisms for recognizing the effort, creativity, 
and achievement of employees. 


e Training. This is listed separately in U.S. EPA (1992), but rigorous training of 
crews who collect data certainly needs to be emphasized in TQM. EMAP makes 
the point that such training should focus on the development of organizational 
and interpersonal skills as well as technical proficiency. 


e Development of a clearly written code of practice for data collection. This 
too is not listed separately under TQM, but it seems to fit there. This code of 
practice is essential for training and for ready reference during data collection. 


Methods 


Approach Adopted 
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Specific quality assurance objectives for Forest Health Monitoring are (U.S. EPA 
1992): 


Compatibility of data. For purposes in Region 6, this means compatibility of 
data within and among Forests and with cooperators such as FIA and BLM. 


Meeting specified levels of uncertainty that the users are willing to accept. 
Such statements should be definitive and either quantitative or qualitative. 


Documentation of data collection methods. Such documentation ensures 
that information on data quality, statistical design, algorithms, protocols, and 
analytical procedures are available for comment. 


Verification and validation of data. A systematic approach to verification and 
validation ensures that all data are subjected to basic standards of accuracy 
to verify their authenticity and reliability. 


To facilitate data verification and validation, U.S. EPA (1992) recommends use of 
portable data recorders, electronic data verification, and remeasurement of some 
field plots as follows: 


Crews should remeasure a representative sample of their own plots so that 
estimates of their “within-field” precision can be determined. 


All field crews should measure a series of reference plots so that measures 
of accuracy and “between-crew” precision can be determined. 


Remeasurement, by a check crew, of a representative sample of plots com- 
pleted by field crews. This, too, can measure crew accuracy and between-crew 
precision. 


Based on the specified sampling objectives, design criteria, directions from the 

FS Washington Office, literature, personal experience, and the concern in Region 6 
for a single inventory database with consistent resource estimates, the approach 
summarized below was adopted. 


The general approach to IM adopted by Region 6 is explained in the following series 
of statements. These statements document the principles and key features of the IM 
system. 


Continue the development of Regional Geographic Information Systems (GIS) 
databases for historical vegetation, current vegetation, and potential natural 
vegetation across all FS lands, using all levels of technology including state- 
of-the-art modeling and satellite imagery. 


Pursue developing techniques for continually updating GIS databases. 


Establish a 5.47-kilometer (3.4-mile) square grid across all NFS land regardless 
of administrative classification or forest allocation. This grid is compatible with 
that used by the Pacific Northwest Research Station FIA unit and is expected 

to be adequate to meet FIA sampling error standards for each National Forest in 
Region 6. This grid system was implemented in 1993 on all NFS land. Ecological 
land units (Bailey 1994, USDA Forest Service 1993) should be identified and 
mapped on the sample units. Establishing this grid in the field is expected to 

be completed by the end of 1996. 


Mapping 


e Through partnership with FIA, develop standards and procedures for measuring 
ground sample units to include data necessary to meet information needs at all 
levels of the FS. Coordinate operational inventory activities with FIA, including 
data collection, quality control, data compilation and analysis, and reporting. The 
work for any activity may be accomplished by one or both administrative units or 
by contractors, according to the intra-Agency agreement. 


e For land management planning, Forest plan implementation, and monitoring 
changes in land and vegetation for all NFS lands excluding wilderness, use a 
more intensive grid (four times more intense) of a 2.74-kilometer (1.7-mile) 
square. Using this approach, Region 6 personnel have the option to repeatedly 
poststratify the databases as necessary to provide critical information to forest 
managers. Information will not be constrained by strata definitions adopted in a 
typical prestratified inventory. 


e lf the standard grid (2.74-km [1.7-mi] square) for data acquisition does not meet 
all the vegetation or modeling needs of the forest planning effort (it is a given that 
it will not), install a supplemental 1.37-kilometer (0.85-mile) grid. It is at this stage 
that other data sources may be incorporated into the information base. For 
example, the occurrence of Pacific yew was a specific need dictating a particular 
GIS mapping overlay. Such a grid “intensification” can be tied to any use or map 
as a basis for supplemental sampling. 


e Sample all vegetation on all NFS land. In addition, data needed to estimate 
biomass and carbon loading should be collected. 


e Maintain flexibility in design to facilitate addressing future questions. This is 
accommodated by establishing a self-weighting, equal probability, systematic 
grid of sampling units on the ground. 


e Maintain the ability to formally incorporate existing information of acceptable 
quality from other survey sources, such as special range, wildlife, recreation, 
and soil surveys. 


e Provide a formal framework for designing and implementing special studies. 


e Ensure that (1) the sampling units do not change over time because of sub- 
jectivity in the way they are defined, (2) estimates are statistically valid and 
defensible, and (3) statistical estimates of reliability are provided with the most 
important estimates. 


Use of maps for estimating area—Historically, Region 6 has used several methods 
to derive area estimates for a variety of purposes. These methods, regardless of how 
the data were acquired, ranged from plot expansion, in sampling, to determining area 
from maps. No effort was made to estimate the error associated with any of the area 
estimates or the effect of that error on a particular forest resource estimate. This has 
been a continual concern for forest planners and should be of concern to all forest 
managers. 
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Sampling Frame and 
Plot Design 


12 


Technological changes in the management of natural resources and expansion of 
issues encompassing the total landscape have necessitated that the design of inven- 
tories be flexible enough to be compatible with maps created from a wide range of 
sources; for example, Landsat imagery and aerial photography used for ecological 
units, predictive modeling, and even arbitrary delineations of the land base. This 
flexibility does not make it easier to provide resource estimates or make planning 
and decisionmaking easier or more certain. 


Mapping strategy selected—Various maps of the resources are acceptable, but 
any map used must be evaluated, quantitatively, for accuracy. Methods for assessing 
accuracy of maps are currently being developed (Czaplewski 1992, 1994) and will 

be incorporated into the IM system when available. 


To select a probabilistic sample, several sampling frames and designs are possible. 
All those considered preceded the Robertson (1993) letter and are included here 
for completeness. We first discuss alternative sampling frames and then alternative 
designs for the ground plot installations. 


Sampling frames—tThe three frames considered were: 
e Existing vegetation polygons (similar to historical stand maps). 
e Ecological polygons (also called ecological units [USDA Forest Service 1993}). 


e Systematic grid of primary sampling units (PSUs). Each PSU contains a cluster 
of plots and line transects to subsample the PSU. 


The advantages and disadvantages of each of these frames are given below. 
Existing vegetation polygons (VP)— 

Advantages: 

e Allow easy inclusion of information from other (local) management surveys. 


e Facilitate impact analysis because they comprise units that are more meaningful 
to resource specialists. 


e Provide estimates of variation among polygons, within given classes, and within 
polygons. 


e Provide a basis for grouping similar vegetation polygons into meaningful strata to 
improve efficiency. 


Disadvantages: 


e Likely to be unstable over time, particularly the boundaries. Clearly, vegetation 
changes over time, and there is concern that the vegetative classification system 
may not be repeatable given new remote sensing information (for example, new 
Landsat data) and newly developed models. 


e May be relatively inefficient, from a sampling perspective, because a VP should 
be relatively uniform or homogeneous within itself. Subplots within the VP would 
not sample much variation if the classification is done accurately. 


Ecological polygons (EP)— 
Advantage: 


e Likely to provide a better base for “ecologically based” management because 
they are more consistent with the management planning process, both in 
developing and implementing the plan. 


Disadvantages: 


e Likely to be a lack of objectivity in defining the EPs, particularly in establishing 
their boundaries. This will lead to variability, among field crews and even within 
individuals, in defining boundaries. 


e Possible major change in definitions of EPs over time. This could affect 
boundaries on the ground and pose problems when estimating trends. 


Systematic grid of primary sampling units (PSUs)— 
Advantages: 

e Totally objective. 

e Consistent over time. 


e Likely to be relatively efficient from a sampling perspective; that is, a cluster of 
sampling locations (over a reasonably large PSU) should cover considerable 
variation. 


e Simple in concept and implementation. 

e Provide self-weighting PSUs. 

e Accommodate any form of poststratification. 
Disadvantages: 

e Less compatible with “ecologically based” management. 
e More difficult to accomplish impact analysis. 


e More difficult to include information from special inventories or impact analysis 
surveys because these most likely focus on specific polygons or ecological 
response units. 


e Altering sampling intensities and keeping a consistent grid limits the flexibility 
of the design. 


Plot designs—All five plot designs considered are cluster designs in the sense 
that a larger PSU is defined, at least conceptually, and then is subsampled. The 
size and shape of the several plot designs, as well as the subsampling details, 
differ considerably. The PSUs for these designs are as follows: 


1. One-hectare (2.47-acre) circular plot. 
2. One-hectare (2.47-acre) rectangular plot. 


3. A satellite system of seven subplots in a hexagonal arrangement. 
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4. A satellite system of five subsample points in a symmetric arrangement covering 
2.5 hectares (6.2 acres). 


5. A transect (linear) arrangement of five subsample points with supplemental 
sampling for particular insect, disease, and tree conditions. 


Each of these plot designs is described in detail below, including the proposed plan 
for subsampling the PSU. 


Alternative (1)—Satellite system of five subplots sampling a fixed-area, circular PSU 
of 1 hectare (2.47 acres). 


e Contains five nested, fixed-area, circular subplots with size changing given the 
variable of interest; for example, 


Grass and forbs—3-meter (9.8-ft) radius 


Regeneration—trees <12.7-centimeters (5-in) in diameter at breast height 
(d.b.h.), 3-meter (9.8-ft) radius 


Trees 12.7-25.4 centimeters (5-10 in) in d.b.h., 6-meter (19.7-ft) radius 
Trees 25.4-101.6 centimeters (10-40 in) in d.b.h., 18-meter (59.1-ft) radius 
Trees >101.6 centimeters (>40 in), entire area of 1 hectare (2.47 acres) 

e The subplots are laid out in a symmetrical design. 


e Some response variables, such as tree mortality, can be tallied for the entire 
hectare or any representative subset of the plot. 


e Transects can be laid out from subplot locations; for example, 
5 at about 10-meter (32.8-ft) intervals for down woody material 
5 at about 10-meter (32.8-ft) intervals for shrub cover 


e No substitution or movement of the plot is allowed for changing vegetation 
conditions in the field. 


Alternative (2)—A rectangular, fixed-area PSU subsampled by rectangular subplots 
located linearly within the PSU. 


e PSU has dimensions of 10 by 1000 meters or 20 by 500 meters (32.8 by 
3280.8 ft or 65.6 by 1640.4 ft). 


e Contains fixed-area rectangular subplots with size based on the variable of 
interest; for example, 


Regeneration—trees <12.7 centimeters (5 in) in d.b.h., 5 by 5 meters (16.4 by 


16.4 ft) 

Trees 12.7-25.4 centimeters (5-10 in) in d.b.h., 10 by 10 meters (32.8 by 
32.8 ft) 

Trees 25.4-101.6 centimeters (10-40 in) in d.b.h., 10 by 20 meters (32.8 by 
65.6 ft) 


Trees >101.6 centimeters (40 in) in d.b.h., entire area of 1 hectare (2.47 acres) 


e Five subplots are laid out at regular intervals to collect information on live trees: 
Random starting point 


Plots may be tied to center of rectangle or along one side to allow for 
coincident boundaries (easier to monument) 


e Mortality can be tallied for the entire hectare. 

e Transects can be laid out along the direction of travel 
5 at about 40-meter (131.2-ft) intervals for down woody material 
5 at about 40-meter (131.2-ft) intervals for shrub cover 

e No substitution or movement is allowed for changing conditions. 

e Plot may have a random or systematic orientation. 


e Allows for maximum variability of the response variables correlated with stand 
basal area within the PSU. 


e PSU may end up crossing several vegetation conditions. 
Alternative (3)—A satellite system of seven subplots in a hexagonal arrangement. 
e The center plot controls the plot configuration and substitution rules. 


e Satellite plots are a combination of fixed-area circular and variable-radius plots; 
for example, 


Regeneration—trees <12.7 centimeters (5 in) in d.b.h., 3.3-meter (10.8-ft) 
radius 


Trees 12.7-91.4 centimeters (5-36 in) in d.b.h., 7 BAF (metric; that is, 7 square 
meters/hectare [30.5 square feet/acre]) 


Trees >91.4-centimeter (36 in) in d.b.h., 18-meter (59.1-ft) radius 
e Distances among subplots are predetermined and fixed so route of travel is easy. 
e No fixed area is associated with the plot design. 
e Mortality can be tallied on satellite plots by using same rules as for live trees. 
e Substitution (point rotation) is allowed for changing conditions. 


Alternative (4)—A satellite system of five subsample points, arranged symmetrically, 
covering about 2.5 hectares (6.2 acres). 


e Design is currently operational (Pacific Northwest Research Station-FlA unit in 
California). 


e Subsample plots are a combination of fixed-area and variable-radius plots; for 
example, 


Regeneration—trees <12.7 centimeters (5 in) in d.b.h., 3.3-meter (10.8-ft) 
radius 


Trees 12.7-91.4 centimeters (5-36 in) in d.b.h., 7 BAF (metric; that is, 7 square 
meters/hectare [30.5 square feet/acre]) 
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Trees 91.4 centimeters (36 in) in d.b.h., 18-meter (59.1-ft) radius 


Two line transects per subsample point are included in the plot design for 
measuring residue or physical and chemical properties. Each transect is 
17 meters (55.8 ft) long. 


Mortality is tallied on subsample points under the same rules as for live trees. 


Lesser vegetation is tallied on a 5-meter (16.4-ft) radius circular plot at each 
subsample point. 


Root disease is mapped on a 17-meter (55.8-ft) radius circular plot at each 
subsample point. 


No substitution is allowed for changing conditions. 


Alternative (5)—A transect (linear) arrangement of five subsample points with 
supplemental sampling for particular insect, disease, and tree conditions. 


Subsample points are a combination of fixed-radius and variable-radius plots; for 
example, 


Regeneration—trees <12.7 centimeters (5 in) in d.b.h., 3.3-meter (10.8-ft) 
radius 


Trees 12.7-101.6 centimeter (5-40 in) in d.b.h., 7 BAF (metric; that is, 
7 square meters/hectare [380.5 square feet/acre]) 


Trees >101.6 centimeter (>40 in) in d.b.h., 18-meter (59.1-ft) radius 


Area sampled is about 2 hectares (4.9 acres), although no fixed area is 
associated with the plot design. 


Twenty supplemental points, spaced at regular intervals between the subsample 
points, follow the same sampling rules as for the five subsample points, but only 
for predetermined conditions; for example, root disease, bark beetles, dwarf 
mistletoe. 


Mortality is tallied on subsample points. 
Line transects can be included between points. 


No substitution of subsample points allowed for changing conditions. 


Plot design selected—Beginning in summer 1993, Region 6 installed vegetation 
inventory plots on 5.47- and 2.74-kilometer (3.4- and 1.7-mile) grids across all NFS 
lands in Oregon and Washington. A version of the alternative (1) plot design was 
selected. The plot design is illustrated in figure 1. The main features of this design 
are as follows: 


PSU size (sample unit) is 1 hectare (2.47 acres). 
Each PSU consists of five subplots each representing 1/5 hectare (1/2 acre). 


Each subplot consists of three concentric, fixed-area plots used to limit 
measurement to certain tree diameter classes. 


On each 0.02-hectare (1/20-acre) fixed-area sample at each subplot, an ecology 
plot is installed based on the plant association guides developed for each Forest. 


[133.9 ft 
“ 485.1 ft — 


1-ha (185.1-ft) radius 


0.076-ha (51.1-ft) radius 


Transect (51.1 ft) 0.02-ha (26.3-ft) radius 


0.004-ha (11.8-ft) radius 


Figure 1—Basic ground plot design for one primary sampling unit, Pacific Northwest 
Region Vegetation Inventory and Monitoring System. 


e On each subplot, a planar intersect sampling method is used to measure down 
woody material and ground cover. 


e All plots and modifications of the basic plot design are permanently established 
for the purpose of remeasurement; that is, estimating long-term trends. 


e The locations of all grid intersections in Region 6 of 1.37, 2.74, and 5.47 kilo- 
meters (0.85, 1.7, and 3.4 mi) were generated by using the GIS available in 
the Regional Office and were documented so that the actual locations can 
be repeatedly duplicated if necessary. Every intersection (PSU location) was 
consecutively numbered. This system is based on the original grid developed 
in 1960 by Region 6, Timber Management; the Pacific Northwest Forest and 
Range Experiment Station; and the Division of State and Private Forestry 
(Bernstein 1960). 
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The concentric plot sizes and associated diameter ranges are as follows: 
— 0.004 hectare (1/100 acre) for seedlings <7.4 centimeter (3.0 in) 
— 0.02 hectare (1/20 acre) for trees 7.6 to 32.8 centimeters (3.0 to 12.9 in) 


— 0.08 hectare (1/5 acre) for trees 32.9 to 81.0 centimeters (13.0 to 31.9 in— 
east side) 


— 0.08 hectare (1/5 acre) for trees 32.9 to 121.7 centimeters (13.0 to 47.9 in— 
west side) 


— 1 hectare (2.54 acres) for trees >81.0 centimeters (32.0 in—east side) and 
>121.9 centimeters (48.0 in—west side) and for rare features such as occurrence 
of root rot or tally of mortality 


Some specific attributes of the selected alternative are: 


e This plot design does not allow for rotation or changes in the layout of any 
subplots, which conforms to the mandate of the Chief of the FS (Robertson 
1993). 


e By not allowing for rotation of subplots, any major change in vegetation must 
be mapped as part of the plot information. 


e The plot design is flexible enough to be modified by a project or Forest to 
meet specific data requirements (for example, sampling for mortality, acquiring 
old-growth training sites), and it could have worked well for the Pacific yew 
inventory. The subplot sizes used to furnish the 1993 and 1994 data are being 
examined to determine if some modification is required. However, the 1-hectare 
(2.47-acre) circular plot still will be used as the PSU. 


The following quality assurance procedures will be followed: 


e Aclearly understandable manual will be provided for Region 6 data collection 
crews. 


e Potential customers of the inventory results will be identified and have an 
opportunity to comment on this paper and on proposed and actual products. 


e All FS data collection crews (direct employees) will be given 10 days of training 
in plot establishment, measurement, and inventory procedures. 


e Demonstrations of plot establishment and plot inspection shall be given to all 
contracted data collection crews. 


e All data collection crews must measure successfully two out of three certification 
plots before they are allowed to work on an inventory project. 


e The work of all data collection crews shall be inspected at a rate of one plot in 
five during the entire project and shall be evaluated by the same set of data 
acquisition standards as defined in their contract. 


e All inspected plots shall become a part of the inventory database as well as part 
of the quality control database. 


e Inspection of the work of all field inspectors shall be done through the Regional 
Office with assistance from other Forest inspectors. 


Data Processing and 
Estimation 


Analyses 


Acknowledgment 


References 


Rigorous quality control of the field data will be built into development of the final 
database. The estimation theory given for mapped designs in the “Review of 
Literature” will be used. The necessary computer algorithms are now being 
developed by the U.S. Department of Agriculture, Forest Service, Timber 
Management Service Center, Fort Collins, CO: 


Procedures for analysis will be developed as data from the IM system become 
available. These will be keyed to the available standard statistical packages. 


We acknowledge the contributions of David Born, Research Forester (now retired), 
U.S. Department of Agriculture, Forest Service, Intermountain Research Station, 
324 25th Street, Ogden, UT 84401, an initial member of the committee. 


Bailey, R.G. 1994. Ecoregions of the United States. 2d ed., rev. Washington, DC: 
U.S. Department of Agriculture, Forest Service. 1:7,500,000; Albers equal area 
projection; colored. 


Bechtold, W.A.; Ruark, G.A.; Lloyd, F.T. 1991. Changing stand structure and 
regional growth reductions in Georgia’s natural pine stands. Forest Science. 37: 
703-717. 


Bernstein, D.A. 1960. Using the master sample grid in continuous forest inventory, 
Division of Timber Management. Portland, OR: U.S. Department of Agriculture, 
Forest Service, Pacific Northwest Region. Unpublished report. On file with: Pacific 
Northwest Region, P.O. Box 3623, Portland, OR 97208-3623. 


Birdsey, R.A.; Schreuder, H.T. 1992. An overview of forest inventory and analysis 
estimation procedures in the eastern U.S.—with an emphasis on components of 
change. Tech. Rep. RM-214. [Fort Collins, CO]: U.S. Department of Agriculture, 
Forest Service, Rocky Mountain Forest and Range Experiment Station. 11 p. 


Brewer, K.R.W.; Mellor, R.W. 1973. The effect of sample structure on analytical 
surveys. Australian Journal of Statistics. 15: 145-152. 


Czaplewski, R.L. 1992. Accuracy assessment of remotely sensed classifications with 
multi-phase sampling and the multivariate composite estimator. In: Proceedings 
16th international biometric conference; 7-11 December 1992; Hamilton, New 
Zealand. Hamilton, New Zealand: Ruakura Agricultural Centre, Statistics Section: 
ae. 


Czaplewski, R.L. 1994. Variance approximations for assessments of classification 
accuracy. Res. Pap. RM-316. Fort Collins, CO: U.S. Department of Agriculture, 
Forest Service, Rocky Mountain Forest and Range Experiment Station. 29 p. 


Deming, W.E. 1975. On probability as a basis for action. American Statistician. 29: 
146-152. 


Feinstein, A.R. 1988. Scientific standards in epidemiologic studies of the menace of 
daily life. Science. 242: 1257-1263. 


Forest and Rangeland Renewable Resources Planning Act. Act of Aug. 17, 
1974. 88 Stat. 476, as amended; 16 U.S.C. 1600-1614. 


? Personal communication. 1995. Roy Mita, Forester, U.S. 
Department of Agiculture, Forest Service, 3825 E. Mulberry St., 
Fort Collins, CO 80524. 


19 


20 


oil 


Hahn, J.T.; MacLean, C.D.; Arner, S.L.; Bechtold, W.A. 1995. Procedures to 
handle inventory cluster plots that straddle two or more conditions. Forest 
Science Monograph. 31: 12-25. 


Hall, J.V.; Frayer, W.E.; Wilen, B.O. 1994. Status of Alaska wetlands. Anchorage, 
AK: U.S. Fish and Wildlife Service, Alaska Region. 33 p. | 


Holt, D.; Smith, T.M.F. 1976. The design of surveys for planning purposes. 
Australian Journal of Statistics. 28: 37-44. 


Kendall, P.L.; Lazarsfeld, P.F. 1950. Problems of survey analysis. In: Merton, R.K.; 
Lazarsfeld, P.F., eds. Continuities in social research: studies in the scope and 
methods of the American soldier. Chicago: Free Press. 


Kohl, M. 1994. Statistical design for the second Swiss forest inventory: a concept 
applying aerial photography and terrestrial sampling. Heft 69. Mitteilungen der 
Eidgen. Forschungs austalt fur Wald, Schnee, und Landschaft. Birmensdorf, 
Switzerland: [Publisher unknown]. 141 p. In German. 


Landis, J.R.; Lepkowski, Eklund, J.M.; Stehouwer, S.A. 1982. A statistical 
methodology for analyzing data from a complex survey: the first national health 
and nutrition examination survey. Series 2, 92. [Place of publication unknown]: 
U.S. Department of Health and Human Services, Public Health Service, National 
Center for Health Statistics. 52 p. 


MacDonald, L.E.; Smart, A.W.; Wissmar, R.C. 1991. Monitoring guidelines to 
evaluate effects of forestry activities on streams in the Pacific Northwest and 
Alaska. Seattle, WA: U.S. Environmental Protection Agency, Region 10. 166 p. 


Messer, J.J.; Linthurst, R.A.; Overton, W.S. 1991. An EPA program for monitoring 
ecological status and trends. Environmental Monitoring Assessment. 17: 67-78. 


Mosteller, F.; Tukey, J.W. 1977..Data analysis and regression. Reading, MA: 
Addison-Wesley Publishing Co. 


National Research Council. 1994. Rangeland health—new methods to classify, 
inventory, and monitor rangelands. Washington, DC: National Academy Press. 
180 p. 


Ouyang, Z.; Schreuder, H.T.; Li, J. 1992. A reevaluation of the growth decline 
in Georgia and Georgia-Alabama. Proceedings, 1991 Kansas State University 
Conference on applied statistics in agriculture; 1991 April 28-30; Manhattan, KS. 
[Place of publication unknown]: [publisher unknown]: 54-61. 


Ouyang, Z.; Srivastava, J.N.; Schreuder, H.T. 1993. A general ratio estimator 
and its application to regression model sampling. Annual Institute Statistical 
Mathematics. 45: 113-127. 


Robertson, F.D. 1993. Forest inventory plot rotation. Letter dated June 24, 1993, 
to USDA Forest Service Regional Foresters, Station Directors, and Area Directors. 
On file with: Office of the Chief, U.S. Department of Agriculture, Forest Service, 
P.O. Box 96090, Washington, DC 20090-6090. 


Ruark, G.A.; Thomas, C.E.; Bechtold, W.A.; May, D.M. 1991. Growth reductions 
in naturally regenerated southern pine stands in Alabama and Georgia. Southern 
Journal of Applied Forestry. 15: 73-79. 


Schreuder, H.T.; McClure, J.P. 1991. Modifying forest survey procedures to 
establish cause-effect, should it be done? 10th world forestry congress; 1991 
September; Paris, France. Revue Forestiere Francaise. Hors Série 4, Sect. D.: 
67-78. 


Schreuder, H.T.; Thomas, C.E. 1992. Establishing cause-effect relationships using 
forest survey data. Forest Science. 37: 1497-1525. 


Schreuder, H.T.; Gregoire, T.G.; Wood, G. 1993. Sampling methods for multi- 
resource forest inventory. New York: John Wiley. 446 p. 


Scott, C.T.; Bechtold, W.A. 1995. Techniques and computations for mapping plot 
clusters that straddle stand boundaries. Forest Science Monograph. 31: 46-61. 


Simmons, W.R.; Bean, J.A. 1969. Impact of design and estimation components 
on inference. In: Johnson, N.L.; Smith, H., eds. New developments in survey 
sampling. New York: John Wiley: 601-628. 


Skinner, C.J.; Holt, D.; Smith, T.M.F., eds. 1989. Analysis of complex surveys. 
New York: John Wiley. 309 p. 


Smith, T.M.F. 1983a. A comment: an evaluation of model-dependent and probability- 
sampling inferences in sample surveys. Journal of American Statistical Association. 
78: 801-802. 


Smith, T.M.F. 1983b. On the validity of inferences from non-random samples. Journal 
of Royal Statistical Society Annual. 146: 394-403. 


U.S. Department of Agriculture, Forest Service. 1993. National hierarchical frame- 
work of ecological units. Washington, DC. 14 p. Administrative document. On file 
with: USDA Forest Service, P.O. Box 96090, Washington, DC 20090-6090. 


U.S. Environmental Protection Agency. 1992. Monitoring and research strategy 
for forests—environmental monitoring and assessment program (EMAP). 
EPA/600/A-91/012. Washingon, DC. 


Williams, M.S.; Schreuder, H.T. 1995. Documentation and evaluation of growth 
and other estimators for the fully mapped design used by FIA: a simulation study. 
Forest Science Monograph. 31: 26-45. 


21 


Appendix 1 
Plot Recording Sheet 


22 


IOAOD 
addy 19009 


s¢]Remarks 


XX 


soinjeoy s 
Aye, dnory 
re pug 8e7/%5| 


Damage/ Severity 


uydaq sugeq 
WIPIM UMOID 


58 


Io]oureIq 


a 

A 
3) 
© 
Q 

n 


o0ue\siq <1 


wanunzy|s 


Current 
Vegetation Data 


Tally Guide 


version 1.4 


1 hectare (2.47 ac.) 


=48") West 


Trees (> 


Actual plot recording sheet used to record data in the Pacific Northwest Region 
Vegetation Inventory and Monitoring System. 


-076 ha. (1/5.3 ac.) 
Trees (13-31.9") East 


Stumps (>= 13"dbh) 


-02 ha. (1/20 ac) 


w 


U. 


Ss. 


Stumps (5-12.9") 

Hardwood Clumps 

Indicator Species 
-004 ha. (1/100) 


GOVERNMENT PRINTING OFFICE: 


(1/5.3 ac.; 51.1' r) 
(1/20 ac; 26.3' r) 
(1/100 ac.; 11.8' r) 


Transect 


3" ID + 
*- For multiple damaging agents 


S - Refer to manual for use 


Down Woody 
Cover Class 


1997 - 589-353 / 63005 REGION NO. 


10 


Max, Timothy A.; Schreuder, Hans T.; Hazard, John W.; Oswald, Daniel D.; 
Teply, John; Alegria, Jim. 1996. The Pacific Northwest Region vegetation 
and inventory monitoring system. Res. Pap. PNW-RP-493. Portland, OR: U.S. 
Department of Agriculture, Forest Service, Pacific Northwest Research Station. 
22 p. 


A grid sampling strategy was adopted for broad-scale inventory and monitoring of 
forest and range vegetation on National Forest System lands in the Pacific North- 
west Region, USDA Forest Service. This paper documents the technical details 

of the adopted design and discusses alternative sampling designs that were 
considered. A less technical description of the selected design will be given 
elsewhere. The grid consists of a regular, square spacing with 5.47 kilometers 
(3.4 mi) between grid points. The primary sampling unit (PSU), established at 
each grid sampling point, consists of a circular, 1-hectare (2.47-acre) plot. The 
PSU is subsampled with a set of different-sized fixed-area subplots, as well as line 
transects, to assess all components of vegetation. The design is flexible and can be 
used with many types of maps. The theory of point and change estimation is 
described, as well as estimates of variation that assess the statistical precision of 
estimates. 


Keywords: Sampling, plot design, fixed-area plots, line intersect sampling, monitoring, 
National Forest System, Pacific Northwest. 
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